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Abstract 

We review some techniques employed in the studies of torques due to bodily tides, 
and explain why the MacDonald formula for the tidal torque is valid only in the zeroth 
order of the eccentricity divided by the quality factor, while its time-average is valid in 
the first order. As a result, the formula cannot be used for analysis in higher orders of 
e/Q . This necessitates some corrections in the current theory of tidal despinning and 
libration damping (though the qualitative conclusions of that theory may largely remain 
correct). 

We demonstrate that in the case when the inclinations are small and the phase lags 
of the tidal harmonics are proportional to the frequency, the Darwin-Kaula expansion 
is equivalent to a corrected version of the MacDonald method. The latter method rests 
on the assumption of existence of one total double bulge. The necessary correction to 
MacDonald's approach would be to assert (following Singer 1968) that the phase lag 
of this integral bulge is not constant, but is proportional to the instantaneous synodal 
frequency (which is twice the difference between the evolution rates of the true anomaly 
and the sidereal angle). This equivalence of two descriptions becomes violated by a 
nonlinear dependence of the phase lag upon the tidal frequency. It remains unclear 
whether it is violated at higher inclinations. 

Another goal of our paper is to compare two derivations of a popular formula for 
the tidal despinning rate, and to emphasise that both are strongly limited to the case 
of a vanishing inclination and a certain (sadly, unrealistic) law of frequency-dependence 
of the quality factor Q - the law that follows from the phase lag being proportional to 
frequency. One of the said derivations is based on the MacDonald torque, the other on 
the Darwin torque. Fortunately, the second approach is general enough to accommodate 
both a finite inclination and the actual rheology. 
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We also address the rheological models with the Q factor scaling as the tidal frequency 
to a positive fractional power, and disprove the popular belief that these models introduce 
discontinuities into the equations and thus are unrealistic at low frequencies. Although 
such models indeed make the conventional expressions for the torque diverge at vanishing 
frequencies, the emerging infinities reveal not the impossible nature of one or another 
rheology, but a subtle flaw in the underlying mathematical model of friction. Flawed is 
the common misassumption that damping merely provides phase lags to the terms of the 
Fourier series for the tidal potential. A careful hydrodynamical treatment by Sir George 
Darwin (1879), with viscosity explicitly included, had demonstrated that the magnitudes 
of the terms, too, get changed - a fine detail later neglected as "irrelevant" . Reinstating of 
this detail tames the fake infinities and rehabilitates the "impossible" scaling law (which 
happens to be the actual law the terrestrial planets obey at low frequencies). 

Finally, we explore the limitations of the popular formula interconnecting the quality 
factor and the phase lag. It turns out that, for low values of Q, the quality factor is no 
longer equal to the cotangent of the lag. 

1 Prologue 

When it shall be found that much is omitted, 
let it not be forgotten that much likewise is performed 

Samuel Johnson, 1755 

In his short work "Untersuchung der Frage ... ," known among the historians also as the 
"Spin-Cycle essay," Immanuel Kant (1754) stated that the Moon not only pulls the Earth, 
but also exerts a retarding torque upon its surface; this torque slows down the Earth's rotation 
and lets go only when terrestrial days become as long as lunar months. Although Kant had in 
mind only the ocean tides, not the bodily ones, we may say that, qualitatively, he predicted the 
celebrated 1 : 1 spin-orbit resonance, the pas de deux wherein Pluto and Charon are locked. 

For the first time, the idea of tidal action not being confined only to the fluid portion of the 
planet but affecting also the solid, so as to induce a state of varying strain, was put forward 
by John Herschel (son of astronomer William Herschel), as a minor aside in a paper devoted 
to volcanism and earthquakes (Herschel 1863). The earliest mathematical description of land 
tides in their dynamics was offered by George Darwin (son of naturalist Charles Darwin and 
great-grandson of poet and philosopher Erasmus Darwin). 

Following his predecessors Roche (1849) and Thompson (1863), who had calculated the 
figure of a static tide, Darwin (1879) assumed the Earth homogeneous and consisting of an 
incompressible fluid. To account for dynamics, he also assumed that the viscosity was the 
sole source of the tidal friction. Relying on this model, Darwin (1880, 1908) derived a tide- 
generated disturbing potential expanded into a Fourier series. Substitution thereof into the 
Lagrange-type planetary equations led him to expressions for the time derivatives of the orbital 
elements via partial derivatives of the disturbing potential with respect to the elements. 

An impressive generalisation of Darwin's work by Kaula (1964), and the subsequent flow of 
new concepts and applications (MacDonald 1964; Goldreich 1966a,b; Goldreich & Peale 1966; 
Singer 1968; Mignard 1979, 1980; Touma & Wisdom 1994; Neron de Surgy & Laskar 1997; 
Krasinsky 2002, 2006; Getino, Escapa & Garcia 2003; Ferraz Mello et al 2008; Efroimsky 2008) 
made bodily tides a rapidly developing area of the planetary astronomy. The vast and growing 
volume of the relevant material leaves one no chance to glean it all in one review. Therefore 
we shall concentrate on one special aspect of this research, the tidal torques emerging from the 
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bodily tides. Moreover, we shall dwell solely on the techniques, not applications. 

Although our review will at times be very critical, it should from the beginning be agreed 
that our criticisms are intended in the spirit of the above quotation from Samuel Johnson. 

Along with reviewing the current state of the field, we shall provide some new results of 
our own. Specifically, we shall address the rheological models with the Q factor scaling as the 
tidal frequency to a positive fractional exponential. We shall demonstrate that, contrary to 
the common opinion, such rhcologics do not cause infinities in the expression for the torque. 
We shall also derive an expression for the tidal torque decelerating a terrestrial planet obeying 
such a rheology. (That the reahstic terrestrial bodies indeed obey this class of rheologies has 
been explained in Efroimsky & Lainey 2007.) 



2 Trivia 

In this section, we shall briefly recall how a satellite-generated potential in a point on or inside 
the planet is expressed through the latitude, longitude, and the radial distance of the point. 

Let us begin from the first principles. The dynamics of point masses located at inertial- 
frame-related positions , 

rui p- ^ mi22 -f ' = Pj ~ Pi ' = 1,...,A^ , (1) 

may be conveniently reformulated in terms of the relative-to-the-primary locations 

Ti = foi = Pi - Po , (2) 
Po standing for the position of the primary. The difference between 

«-Eg^ + g!!^ (3) 

and 



amounts to: 



"^i '^ij V" ^ n^j rj ^ {rrii + mo) n _ dUi 

Ui being the potential: 



y< = - + + , (6) 

n 

with the disturbance 

singled out. This disturbing potential acting on mass rrii is generated by the masses rrij 
other than rrii or the primary. It deviates from the Newtonian one by the amendment 
G rrij rj^ fi ■ fj emerging in the noninertial frame associated with the primary. 
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In the simplest case of one secondary, a satellite of mass mi = M*^^ , located at a planeto- 
centric position Vi = r*, will be creating at some point r2 = -R a perturbing potential 



W{R, n = - G Ml, - > ' 

expandable over the Legendre polynomials (for i? < r* ) by means of the formulae 

1 1 °° / \ ' 

^-;^-;.E(;.) ^*o..) (9) 

and 



R ■ r* Rr* cos 7 1 it! „ , , , , 
- ^ = - -A(cos7) , (10) 



I ^ I 3 "J y» 1^ " 



7 being the angular separation between R and r * , subtended at the point of origin, which 
we shall naturally choose to coincide with the planet's centre of mass. Together, the former 
and the latter formulae yield: 

00 

W(R,r-) = -^±(IL) P,(c„s7) , (11) 



1=2 



where we have neglected the / = term — GM*^^,/ r * , because it bears no dependence upon 
R, and in practical problems is attributed to the principal part of the potential, not to the 
one regarded as perturbation. The angle 7 can be expressed via spherical coordinates as: 

R r* 

cos 7 = — = sin0 sin0* + cos0 cos0* cos(A — A*) , (12) 

it r * 

{R, (f) , X) being the planetocentric distance, the latitude, and the longitude of the point 
where the disturbance is experienced; and (r* , 0* , A*) being the spherical coordinates of the 
satellite. It is customary (though not at all obligatory) to reckon the longitudes from a planet- 
fixed meridian, in which case the subsequent formulae for the potential come out written in a 
reference frame co-rotating with the planet. 

A Legendre polynomial of cos 7 , too, can be expressed via the spherical coordinates: 

PKCOS7) = V II 7 ""1; (2 - 5om) Pim{sm(l)) Pim{sm(l)*) cosm(A - A*) , (13) 

substitution whereof into (11) results in 

00 ' I 

W{R, n = - (^) J2 77X^(2 - Som)Pim{sinmm{sin<l>*) cosm(A - A*) . (14) 

1=2 ^ ^ m=0 ^ 

Evidently, this formalism will stay unaltered, if the role of the tide-raising satellite is played by 
the Sun, or by another satellite, or by another planet. (In this case, what we call Mgat will, in 
fact, denote the mass of the Sun, or of the other satellite, or of the other planet.) Likewise, the 
formalism may in its entirety be applied to a satellite regarded as a tidally-disturbed primary, 
the planet being treated as a tide-raising body (and M*^, now standing for the planetary 
mass) . 
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3 The Kaula expansion for a tidal potential 

Kaula (1961) came up with a remarkable formula 



( — j P;(sin0*) [cosmA* + -\/^ sinmA* ] = 



i+i °° oo ( ( * Q*\ ^ I — i" ■ ( * ^ 6ven 



where Fimp{i) are the inclination functions (Gooding and Wagner 2008); Gipg{e) are the 
eccentricity polynomials identical to the Hansen coefficients X^j^J^^^^-^' ^''^'^^^ ; the notation 

is used to avoid confusion with the inclination; and the auxiliary combinations Vi^p^ are 
defined as:^ 

v^^p^ = (/ - 2p)u;* + (/ - 2p + q)M* + mQ* . (16) 

This development enabled Kaula (1961, 1964) to carry out a transformation from the tide- 
raising satellite's spherical coordinates to its orbital elements and the sidereal time 6* . These 
elements (the scmimajor axis a* , eccentricity e* , inclination i* , pcriapsc uj* . ascending node 
Q* , mean anomaly Ai* ) are introduced in a frame that is associated with the equator but is 
not co-rotating with it. In terms of these parameters, 

W{R,n--^^f:(§) E |[^(2- Wi^Usm0)^F,^,(^*) £ G,,,{e*) 

1=2 ^ ^ m=0 ^ p=0 q=-oo 







i-m even 


cosmA < 


cos ^ 






sin j 


i-m odd 



{w-m even 
) i-m odd 



or, after carrying out the multiplication of the sine and cosine functions: 



Q* Z-^ \ a* J ^ (/-|-m)! 

1=2 ^ ^ m=0 ^ ' 



(15) 



I oo ^ >.i-m even 

5o„.) P^sin^) E E ^^^'^(^*) Z Kp.-^(A + r)) . (18) 

p=0 q=-oo ^ ) l-vn odd 



^ This definition agrees with that by Kaula (1961, 1964, 1966), but differs from the one by Lambeck (1980) 
who incorporated —m9* into v*^^^ . 
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4 Physical assumptions involved in Kaula's theory 

If the primary is not a point mass, it becomes distorted by potential W{R , r *) . The distortion 
of shape will, in its turn, generate some extra potential perturbation whose calculation is 
complicated by the tide-raising potential (18) evolving in time and having a rich spectrum of 
frequencies. The response of the primary's shape to each of these is different and depends on 
the properties of the planet's material. This is a situation where the linear approach becomes 
most helpful, when applicable.^ 

The linear theory of bodily tides comprises two independent assertions. One is that the 
energy attenuation rate (E) at each harmonic depends solely on the frequency x '^^ 
amplitude Epeak{x) > is not influenced by the rest of the spectrum. This is written down 
as {E{x)) = -xEpeak{x)/Q{x), which is equivalent to AEcydeix) = -'2^Epeak{x)/Q{x), where 
AEcycieix) is the one-cycle energy loss, and Q{x) is the quahty factor. The other assertion is 
that each stationary tidal change of the potential, Wi , inflicts on the planet's shape a linear 
deformation. Each of these deformations, in their turn, amend the potential of the primary 
with an addition proportional to the Love number ki . As known from the potential theory, an 
addition proportional to P;(cos7) must be decreasing outside the spherical primary as r"*^'"*"^) . 
Hence, were the external potential perturbation W static (or, equivalently, were the response 
of the material instant) , the tidal addition to the planetary potential would have assumed the 
form^ 

l+l 



;— o V ^ / 



^ fR\ gm;^ ( R \ ^ {l-m)\ ,^ 

I OO J- N i-m even 

- Som) Pim{sm<f>) J2 Fimpin J] { vLp, ' m{X + 9*) ) . 

p=0 q=-oo ^ 'l-m odd 

it! being the mean equatorial (equivolumetric) radius of the planet, R — {R, (f), X) being a 
particular surface point, and r — (r , , A) being an exterior point located right above the 
surface point R, at a. planetocentric radius r > R. 

As we intend to study the effect of this potential on another external body, a similar 
transformation should be applied to the coordinates (r , , A) , to express W through the 
orbital elements of this body. Employment of (15), this time not for r* but for r, leads to: 



^ For most materials, departure from linearity becomes considerable when the strains approach 10 ^ . 
(Karato 2007) 

^ Following MacDonald (1964) and Singer (1968), wo denote the tide-raising potential with W and the 
bodily-tide one with U . In his original paper, Kaula (1964) called these potentials U and T , while in the 
book he switched to U and Uj, (Kaula 1968). Be mindful that we are using a sign convention different from 
that of Kaula. As our forces are negative gradients of potentials, our potentials are negative to those of Kaula. 
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(20) 

I OD I OO 

p=0 q=—oo h=0 j=—oo 

a formula that generalises the tidal theory of Darwin (1908, p. 334) to I and \q\ larger than 
2 . Both Kaula (1964), who derived this milestone result, and Darwin, who had developed its 
simplified version, realised that this machinery would work only after the material's delayed 
reaction to perturbation (18) is somehow taken into account. Until then (20) remains idealised, 
in that it corresponds to an unphysical case of instantaneous response. 

To account for damping, Kaula (1964) followed the path of Darwin (1880, 1908): he en- 
dowed each term of the Fourier series (20) with a real phase lag of its own, eimpq , whereafter 
the ultimate form of Kaula's expansion became 

Uirl f ^V"^ (^)± (2- 

(21) 

( OO ( OO 

p=0 q=~OD h=0 j=—oo 

This empirical method of including dissipation into the picture contains in itself an important 
omission, of which Sir George Darwin was aware, but which was overlooked by his successors. 
Briefly speaking, even in a linear system a dissipation process is not fully accounted for by 
amending phases of the Fourier components. This observation happens to be of relevance in 
the theory of tidal torques. We shall return to this point in section 9. 



5 The two sidereal angles 



Kaula's construction contains a seemingly redundant fixture, which turns out to be an im- 
portant and useful acquisition. This is Kaula's introducing two sidereal angles instead of one. 
As these angles, 9 and 9* , are not orbital elements of the tide-raising and tidally disturbed 
moons, but are parameters characterising the instantaneous attitude of the planet, it may look 
strange that Kaula (1964) assumed them to be different entities. To understand his point, let 
us trace the physical origin of the phase lag. The material of the primary is being deformed 
by a tidal stress whose spectrum contains an infinite number of frequencies, the reaction of the 
material to each of these being different. In a linear regime, the strain has the same spectrum, 
with each harmonic delayed by its own time lag Ati^pq ■ Singer (1968), and later Mignard 
(1979, 1980), assumed that all Atimpq are equal to one another: Atimpq = At. If this were true, 
then in Kaula's series each argument vf^^^ — m9* would have to be substituted with 



J, (delayed) 

Impq 



mO 



^ {delayed) 



Impq 



{t - At) - m9*{t - At) = v^^Jt) - m9*{t) 



Impq 



m9* 



At 



— V 



lmpq\ 



(t) - m9*{t) - {l-2p) uj* + {l-2p + q) M* + m {Q* - 9*) At 



(22) 
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In reality, however, the time lag is a function of frequency, for which reason the delays Atimpq 
will be different for each harmonic involved. This is why the arguments ff^p^ — m9* at the 
moment t should rather be replaced with 



^ {delayed) 

hmpq 



m 



(delayed) 

Impq 



(l - 2p) CO* + {l-2p + q) M* + m (Q* - 9*) At 



(23) 



Specifically, 



V 



jj, (delayed) 

Impq 



Impq 



{I - 2p) dj* + {l-2p + q) M* + m fl* 



''Impq 



At 



Impq 



and 



9] 



(delayed) 



9* - 9* At,. 



'Impq — 1/ L^Hmpq , 

9* being the planet spin rate. In brief, (23) can be rewritten as 



■^Impq 



— m9 



Impq 



Hmpq 



^ 9 ^Impq At Impq 



We see that the total phase lags eimpg introduced by Kaula are given by 



^Impq 



(/ - 2p) 1:0* + {l-2p + q)M* + m (Q* - 9*) At 



''Impq 



^lm,pq Atlmpq 



''Impq 



the tidal harmonic uji^^^ being introduced as 



UJ 



Impq 



{I - 2p) u* + {l-2p + q) M* + m (Q* - 9*) 



the positively-defined physical frequency 

X^lmpq^ l^Lpql = \ {I - 2p) + {I - 2p + q) M* + TU {Cl* - 9* 



(25) 



(26) 



being the actual physical Impq tidal frequency excited in the primary's material. The ap- 
propriate positively-defined time delay At impq depends on this physical frequency, for which 
reason the delays Atimpq are, generally, different from one another. ^ 

The sign on the right-hand side of (24) is simply the sign of uj^^p^ . The sign evidently 

depends on whether m9 falls short of or exceeds the linear combination (l — 2p) Co* + {I — 
2p + q) M* + mQ* ^ {l-2p + q) M* . 

The origin and meaning of the phase lag ei^pq being now transparent, one may express the 
cosine functions in (21) either as 



cos [ {v*^pq - m9*) - {vimhj - m9) - eimpq] 
(where 9* and 9 are identical and cancel one another), or simply as 



cos 



^ (delayed) 

hmpq 



m 9 



, (^delayed) 
Impq 



) - {vimhj - m9) 



(27) 



(28) 



In (28) we have the delayed siderial angle, 6'*^^"'""*' ; separated from the actual angle, 9 , by 
- 9 Atimpq , the time lag Atimpg being a function of ximpq = \ ^^impq \ ■ 



(24) 



* When Kaula was developing his theory, the functional form of the dependence Af(x) was not yet known. 
Reliable data became available only in the final quarter of the past century. See formula (83) below. 
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6 The Darwin-Kaula-Goldreich expansion 
for the tidal torque 

Now we are prepared to calculate the planet-perturbing tidal torque. Since in what follows we 
shall dwell on the low-inclination case, it will be sufficient to derive the torque's component 
orthogonal to the planetary equator: 

r = - M„, ^ , (29) 

Msat being the mass of the tide-disturbed satellite, and the "minus" sign emerging due to 
our choice not of the astronomical but of the physical sign convention. Adoption of the latter 
convention implies the emergence of a "minus" sign in the expression for the potential of a 
point mass: — GM / r . This "minus" sign then shows up on the right-hand sides of (6 - 8) 
and, later, of (19 - 21). It is then compensated by the "minus" sign standing in (29). 

The right way of calculating dU{r)/d6 is to take the derivative of (28) with respect to 6 , 
then to insert (23) into the result, and finally to get rid of the sidereal angle completely, by 
imposing the constraint 9* — 9 . This will yield:^ 

1=2 ^ ^ ^ m=0 ^ ' p=0 q=-oo 

I oo 

^Fi^h{i) Gihj{e) sm[v;^pg - vimhj - Qmpg] , (30) 

h=0 j=—OD 

In the case of the tide-raising satellite coinciding with the tide-perturbed one, Mgat = ^tat i 
and all the elements become identical to their counterparts with an asterisk. For a primary 
body not in a tidal lock with its satellite,^ the torque (30) can be split into two parts. The 
first part is constituted by those terms of (30), in which indices (p, q) coincide with (/i, j) , 
and therefore all vimhj cancel with v^^^^ , provided the tidally-perturbed satellite and the 
tide-raising one are the same l)ody. This component of the torque is, therefore, constant. The 



^ Formally, one can as well differentiate (27) instead of (28), first ignoring the fact that 9* and 9 are 
identical and then, after differentiation, permitting them to cancel one another. Though this method produces 
the same result as the rigorous calculation, it nonetheless remains a formal procedure lacking physics in it. 

^ This caveat is relevant, because in resonances expressions (33 - 35) will require modifications. For example, 
the sidereal angle of a satellite tidally locked in a 1:1 resonance will be: 9 = ^ + uj-\-M.-\- 180° -\-a-\- 0{i'^) , 
letter a denoting the librating angle, which is subject to damping and therefore is normally small (less than 
2" for the Moon). Inserting the said formula for 9 into the expression (25) for the tidal harmonic, we obtain, 
in neglect of —ma : 

^impq = {l-^p-m)^;* + {l-2p + q-m)M* . 

We now see that, since ^ is a function of the other angles, different sets of the indiccs's values will correspond to 
one value of the tidal frequency. We shall illustrate this by considering the so-called anomalistic modes ±A1 in 
the potential. These modes, corresponding to the physical frequency \M \ , are given by [Impq) ~ (201, ±1) 
and also by {Impq) = (220, ±1). Although the m = terms enter the potential, they will not be in the 
torque, as can be observed by differentiating (19) with respect to A, or by differentiating equation (21) with 
respect to —9 , or simply by noticing the presence of the factor m on the right-hand side of (30). Nonetheless, 
we see that there exists a pair of m = 2 terms, which provides an anomalistic input into the torque. This way, 
the case of libration deserves a separate consideration, as it is more involved than that of tidal despinning. 
Specifically, in the case of libration a value of the tidal frequency may correspond to different sets of the indices' 
values. 
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rest of the total sum (30) will be denoted with r . It is comprised of the terms, in which the 
pairs (p , q) differ from {h , j) . Accordingly, these terms contain the differences 

* 

''^Impg ~ '^Imhj — 

{I - 2p + q) M* - {I - 2h + j)M + m{n* -Q) + I {uj* - uj) - 2puj* + 2huj . (31) 

When the tidally-perturbed and tide-raising moons are the same body, this becomes 

^Lp, - ''^imhj = {2h - 2p + q - j)M + {2h - 2p)u , (32) 

whence we see that the oscillating component of the torque, f , consists of two parts. The part 
with h — p — and q — j ^ consists solely of short-period terms, and it averages out 
trivially. 

The mixed-period part of (30), with h — p ^ , consists of both a short-period contri- 
bution dependent upon the mean anomaly, and a long-period contribution depending upon 
the argument of pericentre. All such terms contain multipliers like F2mp{i) -^2m?i(0 ' where 
/i ^ p , and F220 = 3 + 0{i^) , F210 = 3/2 sin i + 0{i^) , F211 = -3/2 sinz + 0{z^) , 
-^221 = 3/2 sin^ i , the other relevant -F2mn's being of higher order than 0{P) . So the 
only long-period terms that we have to consider in (30) involve products: -f2io(0 -^2ii(0 ' 
-^211(0 -^210(0 , -^220(0 -^221(0 , and -F22i(0 -^220(0 • However these products are of order 
O(z^) . Thus, while the short-period terms in (30) average out over one rotation period of 
the moon about the planet, the long-period terms arc of order 0{e^i^) , the coming from 
the G2pq functions. (Indeed, when h and p differ by 1 , then q and j must differ by 2 , 
to eliminate the mean anomaly, i.e., to make the term long-period and not short-period.) So 
both the short- and long-period contributions may be neglected in our approximation.^ 

Had we tried to expand our treatment to higher inchnations, our neglect of the short-period terms would 
remain legitimate, for they still would average out over one rotation period of the satellite about its primary. 

As for the long-period terms, it would be tempting to say that these average out over the apsidal-precession 
period. The latter is much shorter than the time scale of the planetary spin deceleration, a circumstance that 
may seem a safe justification for the neglect of the long-period terms also for higher inclinations. However, a 
word of warning would be appropriate here. As well known from Kozai (1959a), who took into accoimt the 
primary's nonsphcricity, the pericentre of a satellite inclined by about 63° or 117° will neither advance nor 
retard, at least within the first-order (in J2 ) perturbation theory. (For a critical review of Kozai's theory see 
Taff 1985.) Kozai's original attempt to introduce corrections owing to J3 and J4 was flawed because in the 
vicinity of the critical inclinations these terms should be considered not as higher-order but rather as leading. 
His later analysis demonstrated that at these inclinations the satellite's perigee should librate about 90° or 
270° (Kozai 1962). Under these circumstances, the long-period terms in our expression for the torque will 
not be averaged out. We however may neglect this possibility, because in the current work we consider only 
low-inclined moons. 

Another situation, which we exclude from our treatment, is libration of the satellite's periapse about 90° 
or 270° , caused by the pull of a third body (the star or some large neighbouring planet). The possibility of 
such librations may be derived from the presence of the cos 2uj term on the right-hand side of the equation 
for dw/dt in the theory of Kozai (1959b, 1962) - for an easy introduction into this theory see Innanen et al. 
(1997), and for its generalisation to finite obliquities see Gurfil et al. (2007). An important special case of the 
theory of the third-body-caused librations is the one of the satellite getting into a resonance with the third 
body. An indication that such resonances may cause the satellite's uj librate comes from the mathematically 
similar theory of Pluto- Neptune resonances (Williams & Benson 1971): being in resonance with Neptune, Pluto 
has its periapse librating due to a high inclination. (To be exact, the behaviour of Pluto's periapse is dictated 
not only by Neptune, but by the combined influence of all of the four gas giant planets. However, this does not 
change the main point: the outer body or bodies can cause apsidal libration.) 

In the Solar system, none of the large satellites is so highly perturbed as to have a periapse librating around 
90° or 270° due to the above two mechanisms. In theory, though, this remains an option for exoplanets. 
Either librating mechanism might apply also to satellites of minor planets. In our current paper we do not 
consider such moons. 
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Thus we arrive at: 

T = J]2 A;, G M,„i J] X, ^ XI ^;mp(0 XI ^IPI^^") ^^^^Irnpq + f , (33) 

1=2 O m=0 ^ p=0 q=-oo 

the sum standing for the constant (A4-independent) part of the torque, and f denoting the 
oscillating part whose time-average is zero. 

As we pointed in the end of section 5, the sign of the phase lag eimpq depends on whether 
m 9 falls short of or exceeds the linear combination {I — 2p) u* + {l — 2p + q) Ai* + m Q* fa 
(/ — 2p + q) M.* . Now we also understand that, outside resonances, the Impq component of 
the tidal torque experienced by the planet is decelerating if the values oi mO exceed the given 
combination, and is accelerating otherwise. 

Expression (33) gets considerably simplified if we restrict ourselves to the case of / = 2 . 
Since < m < / , and since m enters the expansion as a multiplier, we see that only 
m — 1, 2 actually matter. As < p < / , we are left with only six relevant F's, those 
corresponding to (Imp) = (210), (211), (212), (220), (221), and (222). By a direct inspection 
of the table of Fimp we find that five of these six functions happen to be 0{i) or 0{i'^) , the 
sixth one being F220 = f ( 1 + cos ^) = 3 + O(i^) . Thus, in the leading order of i , the 
constant part of the torque reads: 

^,=2 ^ 2 ^ sin6^^^^ +0{^'/Q) . (34) 

q=—oo 

This is what is called Darwin-Kaula-Goldreich torque, or simply Darwin torque. The principal 
term of this series is 



3 ^ , „5 

2200 2 



= n^^satk2R a ^ sin £2200 ■ (35) 



Switching from the lags to quality factors via formula^ 

Qlmpq — I cot Cijnpq \ , (36) 

we obtain: 

sgncj, sgna;, sgna;, 

I I Impq ^ Impq " Impq , z^//^— 3\ 

sme, = sm e, sgnw, = , = , = = — — + 0{Q ^) 

Impq Impq' Impq /l i ^^+2 ^ / 2 O 

^/-^ + ^Impq Jl + Q ^^^P1 



Impq 



whence 



3 5 2 sgna;oonn 

q=-oo ^220q 



^ The phase lag eimpq is introduced in (23 - 24), while the tidal harmonic uiimpq is given by (25). The 
quality factor Qimpq = \ cot eimpq \ is, for physical reasons, positively defined. Hence the multiplier sgnivimpq 
in (37). (As ever, the function sgn(a;) is defined to assume the values +1 , — 1 , or for positive, negative, or 
vanishing x, correspondingly.) 

Mind that no factor of two appears in (36 - 37), because e is a phase lag, not a geometric angle. 
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Now, let lis simplify the sign multiplier. If in expression (25) for cj^^^^^ we get rid of the 

redundant asterisks, replace^ M. with Mq + n ~ n, and set / = m = 2 and p — 0, the outcome 
will be: 



2uj + {2 + q)n + 2n- 29 



= sgn 



u + [1 + 



n + n- 9 



As the node and periapse precessions are slow, the above expression may be simplified to 

q 



sgn 



All in all, the approximation for the constant part of the torque assumes the form: 



1=2 



2 E ^^sat R a-' Clqie) k, sgn 



(l + I) n-9j+ O(zVQ) + 0(Q-3) 



That the sign of the right-hand side in the above formula is correct can be checked through the 
following obvious observation: for a sufficiently high spin rate 9 of the planet, the multiplier 



sgn 



1 + 



n - 9 



becomes negative. Thereby the overall expression for r acquires a 
"minus" sign, so that the torque points out in the direction of rotation opposite to the direction 
of increase of the sidereal angle 9 . This is exactly how it should be, because for a fixed q and 
a sufficiently fast spin the ^'s component of the tidal torque must be decelerating and driving 
the planet to synchronous rotation. 

Expansion (38) was written down for the first time, without proof, by Goldreich & Peale 
(1966). A schematic proof was later offered by Dobrovolskis (2007). 



7 The MacDonald expression for the tidal torque 

The idea of representing the tidal pattern with one bulge is often mis- attributed to MacDonald 
(1964) who in fact borrowed it from Gerstenkorn (1955). This approach was furthered by 
Singer (1968) and greatly advanced by Mignard (1979, 1980). It is the latter two authors who 
realised that Gerstenkorn's single-bulge simplification was acceptable only with a frequency- 
independent At, not with a frequency-independent Q as in MacDonald (1964). Nevertheless 
we shall call this approach "the MacDonald torque" , to comply with the established convention. 

^ While in the undisturbed two-body setting M. = M.o + n{t — tQ) and M. = n, under perturbation these 
relations get altered. One possibility is to introduce (following Tisserand 1893) an osculating mean motion 
n{t) = ^ ii/a{tY , and to stick to this definition under perturbation. Then the mean anomaly will evolve as 

M = M„{t) + J n(t) dt , whence M = M^it) + n{t) . 

Other possibilities include introducing an apparent mean motion, i.e., defining n either as the mean-anomaly 
rate dA4/dt, or as the mean- longitude rate dL/dt = dfl/dt + duj/dt + dM./dt (as was done by Williams et 
al. 2001). It should be mentioned in this regard that, while the first-order perturbations in a{t) and in the 
osculating mean motion ii/a{t)^ do not have constant parts leading to secular rates, the epoch terms typically 
do have secular rates. These considerations explain why there exists a difference between the apparent mean 
motion defined as dL / dt (or as dA4 /dt) and the osculating mean motion ii/a{t)^ . 

In many practical situations, the secular rate in Aio is of the order of the periapse rate, while the secular 
rate in Lq turns out to be smaller. Hence the advantage of defining the apparent n as the mean-longitude rate 
dL/dt, rather than as the mean-anomaly rate dM/dt. (At the same time, for a satellite orbiting an oblate 
planet the secular rates of Mq , Lq , and periapse are of the same order.) 

Although the causes of orbit perturbations are beyond the scope of our paper, we would mention that in the 
expression (25) for uj^ the notations M. ,d) , and O generally imply the secular rate. 
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For the same reason, the afore-described Darwin-Kaula-Goldreich expansion will be referred 
to simply as "the Darwin torque" . 

In the preceding section, the Darwin torque's component orthogonal to the equator was 
conveniently given by the fundamental formula (29). Within the MacDonald approach, it will 
be more practical to write the torque as a derivative taken with respect to the longitude. The 
torque acting on the tidally disturbed satellite of mass Mgat is — Mg^t dU/dX, while the 
torque that this moon exerts on the planet is this expression's negative: 

r(fO = . (39) 

Speaking rigorously, the formula furnishes the torque's component perpendicular to the plan- 
etary equator. As can be seen from (44), formula (39) coincides with (29) for low inclinations. 



7.1 Simplifications available for low i 

In principle, we can as well insert into 



= E -) ^^(^' ^ 

1=2 V ^ / 



) 



the "raw" expression (14), the one as yet "unprocessed" by (15). This will give us 

Uir) = -GM:„,E^^^T^^E77-^(2-5ojP,^(sin0)P,Usin0*) cosm(A-A*) 
1=2 r r* rn=o^''^^^- 

or, for low inclinations of both the tidally-perturbed and tide-raising satellites: 

' (J- V 

U{r) = -GM:,, E i+i i+i E 77tS^2 - Som)PUO)PUO) cosm(A - A*) 
1=2 r r* ^^qK''^^)- 

+ 0(1^) + 0{i*^) + 0(11*) 

At this point we once again are faced with the question of how to bring damping into the 
picture, i.e., how to take care of the delayed reaction of the planet's material to the tidal 
stress. It is tempting to substitute mX with its delayed value. Then instead of cosm(A — A*) 
we get 



(delayed) 

cos I m A — m A 1 = cos ( m X 



mX — m X 



(42) 



This trick, suggested by Kaula (1968, page 201),^" has a physical justification only if At is the 
same for all frequencies, a model pioneered by Singer (1968) and furthered by Mignard (1979, 



Mind the difference in notations. While in the original paper Kaula (1964) denoted the phase lags with 

^impq ' book Kaula (1968) called them </3;„j,^ . For the longitudinal lag 2 X At emerging in our formula 

(46), Kaula (1968) used notation 26 . This way, in the terms used by Kaula (1968) in his book, the geometric 
angle subtended at the primary's centre between the directions to the bulge and the moon is called 2 5 , not 5 
as in most literature. 
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1980). It can be shown that this model is equivalent to the following rheological law:^^ 

1 



Q 



Impq 



(43) 



Even then, though, it remains unclear how to connect the longitude lag m \ At with one 
or another Qi„^pg , in the spirit of (37). To see what can be done, write down the longitude 
(reckoned from a fixed meridian on the rotating planet) as 

X^-0 + Q + u; + u + O(i^) ^-e + n + u; + M + 2e sin Al + O(e^) + 0{i'^) , (4 

p being the true anomaly. Thence, in neglect of the nodal and apsidal precessions, the cosine 
becomes: 



cos 



m X — mX I + m X At 
or, equivalently: 



] + m A *At ^ = cos m A - mX* ] + m z>* - r] At ^ , (45) 



cos 



mX — mX 



+ m X At]= cos 



mX — mX 



+ m 



n* -e* 



At + 2me*n*At cosM* + 0{e^) 



Insertion of (45) into (41), along with substitution of r*{t) by r*{t — At) , leads us to 

(< — my. 



u{r) = -gm;,, 



R 



1=2 r 



{t)'^\*{t-Atf' t:^,il + rn)\ 



(2 - 8om)Plm{Q)Plm{^) COS m 



A- A 



+ m 



0* -e* 



At +o(i2) + o(r^) + o(M*) 



If we take into account only the / = 2 contribution, expression (41) will simplify to 

5 2 

m = - E 1^(2 - 5om)P2m{mm{^) COSm{X - A*) + 0{f) + 0{e') + 0(1^) 

r{t) r [t) rn=o ^ 

where only the m = 2 term is important. In the presence of dissipation, the appropriately 
simplified version of (48) will read: 



4 r{t) r'{t-Aty 



cos 



2 A - 2A 



+ 2 



At) + 0{i'') + Oii*^) + Oiii*) , 



Combining (24) with the relation Q = 1/ tan|e| . we sec that setting all ^t^^^^ equal to the same At is 
equivalent to saying that the quality factor scales as the inverse frequency: Q = l/(x At) , provided, of course, 
that the Q factor is large. As can be seen from (37), a more exact relation will read: sin(x At) = + Q'^ , 

so that X At = Q-^ + 0{Q~^) . 

Very special is the case when the values of the quality factor are very low (say, much less than 10). In this 
situation, the interconnection between the quality factor and the phase lag becomes quite different from the 
customary formula Q = cot e . Sec the Appendix for details. 

In (48), we may neglect the A-independent term with m = , because our eventual intention is to find the 
torque by differentiating U{r) with respect to A. We may also omit the m = 1 term, because P2i(0) = . 
This omission brings up an error of order 0{ii*) into equations (41), (47 - 50), and (56) 
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while the corresponding expression for the torque exerted by the sateUite on the planet will, 
in this approximation, be given by 



r(r) = M 



sat' 



dUjf) ^ 3GM:^,M,atk2R 
9^ 2r{tfr'{t-Atf 



sm 



2A-2A 



+ 2 



At) +0(^') + 0(^*') + 0(^^ 



In the case when the tidally disturbed satellite coincides with the tide-raising one, i.e., when 
A = A* and Mgat = M*^^. , we obtain: 



r = IcM^^h ^ . f ^ . sin (2 
^ r{t) r{t - At) 



u- 9 



At) + O(zVg) 



(51) 



lGM^,,hjSm (2 



0-9 



At)+0{iyQ) + 0{en/Q\) , 



where the error 0{en/Q^x) emerges when we identify the lagging distance r{t — At) with 
r = r{t) . Replacement of r{t — At) with r is convenient, though not necessary. In subsection 
7.2 below, we shall explain that, after averaging over one revolution of the moon about the 
planet, the error caused by this replacement reduces to 0{^t? jQ^^) , which will be less than 
the largest error. 

The MacDonald torque (51) is equivalent to the Darwin torque (34) with an important 
proviso that all time lags At^^^^ are equal to one another or, equivalently, that the rheological 
model (43) is accepted. Physically, the special case of equal time lags is exactly the case when 
the tide may be rigorously interpreted as one double bulge of a variable rate and amplitude. 
Mathematically, this model enables one to wrap up the infinite series (34) into the elegant 
finite form (51). Formally, this wrapping can be described like this: expression (51) mimics 
the principal term of the series (34), provided in this term the multiplier G|qo is replaced with 
unity, a is replaced with r , and the principal phase lag 



S200 = 2 (n - 9) At 

is replaced with the longitudinal lag or, possibly better to say, with the quasi-phase 



9) At 



(52) 



(53) 



Thus we see that within the MacDonald one-variable-bulge formahsm the longitudinal lag 
(53) is acting as an instantaneous phase lag associated with the instantaneous tidal frequency 
X = 2 |z> — ^1 . This is why we may call it simply e , without a subscript. Evidently, e is (up 
to a sign) twice the geometrical angle subtended at the primary's centre between the directions 
to the moon and to the bulge. 

The geometric meaning of the longitudinal lag being clear, let us consider its physical 
meaning, in the sense of this lag's relation to the dissipation rate. For some fixed frequency 
Ximpg ) the corresponding phase lag ei^pg is related to the appropriate quality factor via 



An attempt to generalise this simplified approaeli to arbitrary inclinations was undertaken by Efroimsky 
(2006). While for constant time lags that generalisation is likely to be acceptable, it remains to be explored 
whether it offers a practical approximation for actual rheologies (60). 

As the subtended angle is \{u — 6) At\ , its double is equal to the absolute value of e , and not to that 



of 



'220g 



2{n-9)At. 



15 



^/Qimpq = idJi \ eimpq \ ■ To keep the analogy between the true lags and the instantaneous 
lag (53), one may conveniently define a quantity Q as the inverse of tan |e| . This will enable 
one to express the MacDonald torque as 



5 

3 ^ , ,0 , R 



o GMl, h—-^ ^ sin 6 + 0{{'/Q) 

^ r{t) r{t-At) 

(54) 



5 

3 ^ , ,9 1 R 



= - GMl, h^Q-' sgn(z> -6) + 0(^Vg) + 0{en/Q'x) + 0{Q- 



■3\ 



2- sac 6 



Since Q was defined as 1/ tan|e| , it is not guaranteed to deserve the name of an overall 
quality factor. At each particular frequency Ximpq ' corresponding quality factor Qimpq = 
1/ tan |e/^pg| is related to the peak energy of this mode, Ep^akiXimpq) ' one-cycle 
energy loss at this frequency, ^E^ydeiXimpq) ' via 

A£;<,*(X,„„) = - '"''r'^'"-'"' . (55) 

^ Impq 

However, it is not at all obvious if the quantity Q defined through the longitudinal lag as 
Q = l/tan|e| interconnects the overall tidal energy with the overall one-cycle loss, in a 
manner similar to (55). The literature hitherto has always taken for granted that it does. 
However, the proof (to be presented elsewhere) requires some effort. The proof is based on 
interpreting x = 2 |z> — ^| as an instantaneous tidal frequency. 

The interconnection between Q = 1/ tan|e| and the overall energy-damping rate mimics 
(55) only up to a relative error of order 0{en/Qx) = 0{enAt) , i.e., up to an absolute error 
of order 0{eQ~^n At) . This is acceptable, because in realistic settings nAt -C 1 . 

7.2 Further simplifications available in the zeroth order of en/Qx 

Suppose we ignore the difference between r and r*, which arc the two locations of the same 
satellite, separated by the time lag owing to the tidal response. We shall now demonstrate that, 
though the relative error of this approximations is 0{en/Qx) , after averaging over a satellite 
period this approximation brings only a 0{e^n^/Q'^x^) relative error into the expression for 
the torque. 

From the well-known formulae r = a(l — e^)/(l -|- e cosi/) and dv/dM = (1 -|- 
e cosz/)^/(l — e^)^/^ we see that 

Ar^r{t)-r{t-At) = -0^-^^^ s^n u Au + O {e {Auf ) = - n At + O {e {n Atf ) 

The time lag is interconnected with the phase shift and the quality factor via the relations 

xAt = e ^ , 
X = 2 1^ — z> I being the instantaneous tidal frequency. Hence 

Ar = r(t) — r(t — At) ~ — a — — sin i/ . 

Q X 
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As Ar is proportional to sin i/, only terms quadratic in Ar survive averaging. Thus, while in 

5 

[/ = - ^GM:^,k2 ^ cos(2A - 2A* + e) + 0{en/Qx) + 0{i^) + 0{i*^) + 0{ii*) , (56) 
and 

5 

T = ^GM,at k2 % sin(2A - 2A* + e) + 0{en/Qx) + , (57) 



2 oui, --sat ^6 



the relative error is 0{en/Qx) + 0{i^) , in the averaged expression 
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2 6 

(t) = - ^^^f^ ( ^ «i^O + 0(eV/gV) + 0(^7g) (58a) 



2 o 4 

— ' r-ZW — 



SGM.^^/csi? 1 /--i? sine di/ +0(eV/gV) +0(zVg) . (58b) 



^4 



it is only 0{e^nyQ\^) + ©(z^) . 

In the above expressions, we asserted after the differentiation that M*^, = Mgat and 
A* = A , implying that the tide-generating and tidally-perturbed moons are one and the same 
body. As soon as A is set to be equal to A* , the sine function in (58) becomes sine l/Q ■ 
So, while the relative error in (58) is 0{e^n^/Q^x^) + 0{i^) , the absolute error becomes 
0(e2nVgV) + 0(^VQ) . 

The error 0{e^n'^ /Q^x^) becomes irrelevant for two reasons. First, our substitution of 
sine with tane = 1/Q generates a relative error 0{Q~^) , i.e., an absolute error 0{Q~^) . 
Second, as explained in the end of subsection 7.1, the uncertainties inherent in our definition 
of the overall quality factor Q entail an absolute error 0{en/Q'^x) ■ Each of these two errors 
exceeds 0{e^v? /Q^x^) . We can then write: 

{r)= - '-^^ { 55^^^^ ^ > + 0«-) + 0(.V« + 0(eu/Q\) (59a) 

3GM^^,k2R 1 /"^^ R sgn{e - u) 



/ ^ ^ ^.^^L^ + 0{Q-') + O(zVg) + 0{en/Q\) . (59b) 



We recall that time averages over one revolution of the satellite about the primary are given by 



(...). J- r ... dM = li-^ r ... , , 

2 TT Jo 2 TT Jo (1 + e cos v) 

while the planetocentric distance is r = a (l — e^) / (1 + e cosz/) , with v being the true anomaly. This way, 

r' . , _ (1 - e^)^/^ i?^ . 



( — ^ sme ) = — / — ^ (jiiiu ^ 

^ r-6 ' 2 7r io (1 + e cosi^f 



sme 



(1 - e^)'^' f'^ if ^.^^ di^ ^ R_ 1 R_ 

2tt Jo ^'""^ a2 (1 - e'^f (1 _ e^)^^^ Jo 



sine du 
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8 Use and abuse of approximation (56 - 59) 



Just as with the formula (56) for the potential, the elegant expression (57) for the torque 
remain correct only to the zeroth order in e/Q, while (59) is vahd to the first order. This 
is the reason why the convenience of this approximation and of its coroUaria is somewhat 
deceptive. Nevertheless, the (56) - (59) were employed by many an author. 

Goldreich & Peale (1966) used them to build a theory containing terms up to . Now 
we see that some coefficients in their theory of capture into resonances must be reconsidered. 
The same pertains to some coefficients in the theory of Mercury's rotation, recently offered 
by Peale (2005). Fortunately, the key conclusions of Peale (2005) stay unaltered, despite the 
corrections needed in the said coefficients.^^ 

Interestingly, Kaula (1968) fell into this temptation, and so did Goldreich (1966b). Equation 
(4.5.29) in Kaula (1968), as well as equation (15) in Goldreich (1966b), is but the above formula 
(59) with the inverse quality factor taken out of the integral: 

Kaula ^ _ ^GM k^R 1^ r " R ^ ^ _ i,) . (Kaula 1968, eqn4.^ 

Besides the afore-mentioned fact that this approach contains a relative error 0{en/Qx} + 

0{Q~^) + 0{i'^) , it suffers a greater defect. Taking Q out of the integral is illegitimate, 
because it implies frequency-independence of Q . This is then incompatible with Kaula's 
implicit assumption of a constant At, an assumption tacitly present in (42).^'^ 

Goldreich (1966b) and Kaula (1968) used this oversimplified formula to investigate libra- 
tions of a satellite trapped in a 1:1 resonance. Other authors used it to evaluate despinning 
rates of bodies outside this resonance. We shall dwell on the latter case in section 10. 



9 Can the quality factor scale as a positive power of the 
tidal frequency? 

As of now, the functional form of the dependence Q(x) for Jovian planets remains unknown. 
For terrestrial planets, the model Q ~ 1/x is definitely incompatible with the geophysical 
data. A convincing volume of measurements firmly witnesses that Q of the mantle scales as 
the tidal frequency to a positive fractional power: 

Q ^ [E^f , where a = 0.3 ± 0.1 , (60) 

E being an integral rheological parameter with dimensions of time. This rheology is incom- 
patible with the postulate of frequency-independent time-delay. Therefore, insertion of the 
realistic model (60) into the formula presented in section 8 will remain insufficient. An honest 
calculation should be based on averaging the Darwin-Kaula-Goldreich formula (38), with the 
actual scaling law (60) inserted therein, and with the appropriate dependence ^tij^pqiXimpq ) 
taken into account (see formula (83) below). 

Stan Pcalc, private communication, 2007. 

This would be incompatible with the MacDonald (1964) treatment as well, because MacDonald's formalism 
necessitates the rheology Q ^ 1/x • We shall return to this point in subsection 10.1. 



18 



9.1 The "paradox" 

Although among geophysicists the scahng law (60) has long become common knowledge, in 
the astronomical community it is often met with prejudice. The prejudice stems from the fact 
that, in the expression for the torque, Q stands in the denominator: 

r ~ 1 . (61) 

At the instant of crossing the synchronous orbit, the principal tidal frequency X2200 b®*^^"^®^ 
nil, for which reason insertion of 

Q ~ , a > (62) 
into (61) seems to entail an infinitely large torque at the instant of crossing: 

J — — >-oo, forv— >-0, (63) 

a clearly unphysical result. 

Another, very similar objection to (60) originates from the fact that the quality factor is 
inversely proportional to the phase shift: Q ~ 1/e. As the shift (24) vanishes on crossing 
the synchronous orbit, one may think that the value of the quality factor must, effectively, 
approach infinity. On the other hand, the principal tidal frequency vanishes on crossing the 
synchronous orbit, for which reason (60) makes the quality factor vanish. Thus we come to a 
contradiction. 

For these reasons, the long-entrenched opinion is that these models introduce discontinuities 
into the equations and can thus be considered as unrealistic approximations for rotating bodies. 

It is indeed true that, while law (60) works over scales shorter than the Maxwell time (about 
10^ yr for most minerals), it remains subject to discussion in regard to longer timescales. 
Nonetheless, it should be clearly emphasised that the infinities emerging at the synchronous- 
orbit crossing can in no way disprove any kind of rheological model. They can only disprove 
the flawed mathematics whence they provene. 

9.2 A case for reasonable doubt 

To evaluate the physical merit of the alleged infinite-torque "paradox" , recall the definition of 
the quality factor. As part and parcel of the linearity approximation, the overall damping inside 
a body is expanded in a sum of attenuation rates corresponding to each periodic disturbance: 

(£;) = ;^(£;(x,)) (64) 

i 

where, at each frequency Xi > 

mx.))^ -2..<||}1.-X.^ . (65) 

( . . . ) designating an average over a flexure cycle, E{Xi) denoting the energy of deformation 
at the frequency Xi > ^-^id QiXi) being the quality factor of the medium at this frequency. 

This definition by itself leaves enough room for doubt in the above "paradox" . As can be 
seen from (65), the dissipation rate is proportional not to 1/Q{x) but to xlQix) ■ This way. 
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for the dependence Q ~ ) the dissipation rate {E) will behave as x^^" . In the limit of 
X ^ , this scaling law portends no visible difficulties, at least for the values of a up to unity. 
While raising a above unity may indeed be problematic, there seem to be no fundamental 
obstacle to having materials with positive a taking values up to unity. So far, such values of 
a have caused no paradoxes, and there seems to be no reason for any infinities to show up. 



9.3 The phase shift and the quality factor 

As another preparatory step, we recall that, rigorously speaking, the torque is proportional 
not to the phase shift e itself but to sin e . Prom (37) and (60) we obtain: 

sin e I = — 1^=^= = — ^^^^=^^^= . (66) 

+ Q \Ji + s X 

We see that only for large values of Q one can approximate | sine | with 1/Q (crossing of 
the synchronous orbit not being the case). Generally, in any expression for the torque, the 
factor 1/Q must always be replaced with 1/ >/l + . Thus instead of (61) we must write: 

11 

T ~ sin e = — , = — , , (67) 

1 + Q + E X 

E being a dimensional constant from (60). 

Although this immediately spares us from the fake infinities at x ^ 0, we still are facing 
a strange situation: it follows from (66) that, for a positive a and vanishing Xi the phase 
lag e must be approaching 7r/2, thereby infiating the torque to its maximal value (while on 
physical grounds the torque should vanish for zero x)- Evidently, some important details are 
still missing from the picture. 



9.4 The stone rejected by the builders 

To find the missing link, recall that Kaula (1964) described tidal damping by employing the 
method suggested by Darwin (1880): he accounted for attenuation by merely adding a phase 
shift to every harmonic involved - an empirical approach intended to make up for the lack of 
a consistent hydrodynamical treatment with viscosity included. It should be said, however, 
that prior to the work of 1880 Darwin had published a less known article (Darwin 1879), in 
which he attempted to construct a self-consistent theory, one based on the viscosity factor of 
the mantle, and not on empirical phase shifts inserted by hand. Darwin's conclusions of 1879 
were summarised and explained in a more general mathematical setting by Alexander (1973). 

The pivotal result of the self-consistent hydrodynamical study is the following. When a 
variation of the potential of a tidally distorted planet, U{r) , is expanded over the Legendre 
functions P;r„(sin0) , each term of this expansion will acquire not only a phase lag but also 
a factor describing a change in amplitude. This forgotten factor, derived by Darwin (1879), 
is nothing else but cose . Its emergence should in no way be surprising if we recall that the 
damped, forced harmonic oscillator 

evolves as 

(-7 + i V''^2-7^ ) * (-7-i-y/a;2_^2 ) f F COS 6 j(At-e) 

xit) = Ci e + C2 e + e 

ujI - y 
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(68) 



(69) 



where the phase lag is 



tan 6 = , (70) 

and the first two terms in (69) are damped away in time.^^ 

In the works by Darwin's successors, the allegedly irrelevant factor of cos e fell through the 
cracks, because the lag was always asserted to be small. In reality, though, each term in the 
Fourier expansions (21), (30 - 35), and (38) should be amended with cose^^^^ . Likewise, the 

correct versions of (50 - 51) and (54) should contain an extra factor of cose^^^Q ■ For the same 
reason, instead of (67), we should write down: 



a a 



T ~ I cose sme| = ; / = ,2a 2a ■ (71) 



At this point, it would be tempting to conclude that, since (71) vanishes in the limit of x — > , 
for any sign of a , then no paradoxes happen on the satellite's crossing the synchronous orbit. 
Sadly, this straightforward logic would be too simplistic. 

In fact, prior to saying that cos e sin e — > 0, we must take into consideration one more 
subtlety missed so far. As demonstrated in the Appendix, taking the limit of Q — > is a 
nontrivial procedure, because at small values of Q the interconnection between the lag and 
the Q factor becomes very different from the conventional Q — cot |e|. A laborious calculation 
shows that, for Q < 1 — 7r/4 , the relation becomes: 



sine cose | = {3Q) 



1/3 



l_4(3g)2/3 + 0(g4/3) 

5 



which indeed vanishes for Q — > 0. Both e22oo and the appropriate component of the torque 
change their sign on the satellite crossing the synchronous orbit. 

So the main conclusion remains in force: nothing wrong happens on crossing the syn- 
chronous orbit, Q.E.D. 



10 Tidal despinning. 

The following formula for the average deceleration rate ^ of a planet due to a tide-raising 
satellite has often appeared in the literature: 

{e ) = - /C \e A{e) - n N{e) ] , (72) 



As demonstrated by Alexander (1973), this example indeed has relevance to the hydro dynamical theory 

of Darwin, and is not a mere illustration. Alexander (1973) also explained that the emergence of the cose 
factor is generic. (Darwin (1879) had obtained it in the simple case oi I = 2 and for a special value of the 
Love number: fcg = 1.5.) 

A further investigation of this issue was undertaken in a comprehensive work by Churkin (1998), which 
unfortunately has never been published in English because of a tragic death of its Author. In this preprint, 
Churkin explored the frequency-dependence of both the Love number ^2 and the quality factor within a broad 
variety of rhcological models, including those of Maxwell and Voight. It follows from Churkin's formulae that 
within the Voigt model the dynamical k2 relates to the static one as cos e . In the Maxwell and other models, 
the ratio approaches cose in the low- frequency limit. 
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where 



^(e) = (l + 3e^ + ^e^) {l - e^Y'" , (73) 

and 

A^(e)= (l + ye^ + f e^+^e«) , (74) 

9 being the sidereal angle, 9 being the primary's spin rate, JC being some constant, and the 
angular brackets designating an average over one revolution of the secondary about the primary. 
This expression was derived by different methods in Goldreich & Peale (1966) and Hut (1981), 
and was later employed by Dobrovolskis (1995, 2007) and Correia & Laskar (2004)^^. 

10.1 Derivation by means of the MacDonald torque 

The following proof of (72 - 74) is implied in Goldreich & Peale (1966) and is presented in more 
detail in Dobrovolskis (2007). Their starting point was the MacDonald torque (59). Hut (1981), 
who approached the issue in the language of the Lagrange-type planetary equations, took into 
account, in the disturbing function, only the leading term of series (21), and substituted the 
principal tidal frequency X2200 ~ 2\6 — n\ with the synodal frequency x = 2\6 — ul . 
Thereby, his approach was equivalent to that of Goldreich & Peale (1966) and Dobrovolskis 
(2007). 

Although not necessarily assumed by these authors, ^° their method, as we saw in the sec- 
tion above, inherently implied the following assertions: 

(I) The quantity % = 2 |6' — z>| is treated as an instantaneous tidal frequency. Accord- 
ingly, the overall quality factor Q is implied to be a function not of the principal frequency 
X2200 instantaneous frequency x ■ 

(II) The functional form of this dependence is chosen as Q — (Ai)~^ > where At is 
the time lag. 

(III) The time lag At is frequency-independent. This assertion is equivalent to (II), as 
can be demonstrated from (24). 

Beside this, those authors neglected the ordcr-en/Qx difference between r and r* in (54), 
generating a relative error in r of order 0{en/Qx) (which, luckily, reduced to 0{e^v? /Q'^x^) 
after orbital averaging) . They also substituted sin e with l/Q , causing a relative error of 

Our formula (78) differs from formula (4) in Correia and Laskar (2004) by a factor of n/x, because in 
Ihid. the quality factor had been introduced as l/(n At) and not as 1/ (x At) . 

It should be mentioned that the original treatment by MacDonald (1964) is inherently contradictory. On 
the one hand, MacDonald postulates (following Gerstenkorn 1958) that there exists one overall double bulge. 
As explained in subsection 7.1 above, this assertion unavoidably implies constancy of the time lag At, so 
that Q ^ 1/x and e ~ x- However, MacDonald (1964) erroneously set Q (and, thence, also e) frequency- 
independent, an assertion incompatible with his and Gerstenkorn's postulate of existence of an overall double 
bulge. 

Whenever in the current paper we refer to MacDonald's torque, we always imply his postulate that one 
double bulge exists. At the same time, to make the MacDonald-Gerstenkorn treatment consistent, we always 
adjust the MacDonald-Gerstenkorn treatment by letting Q and e scale as 1 /x and x ) correspondingly. 



22 



order 0(1/(5^) , because in reality Q is the reciprocal of tane, not of sine. 



Assertion (II) can be written down in more generic notation: 

Q = [Ex)" , with a = - 1 . (75) 

This form of the scaling law is more convenient, for it leaves one an opportunity to switch to 
different values of a. For any value of a (not only for —1), the constant S is an integral 
rheological parameter (with the dimension of time), whose physical meaning is explained in 
Efroimsky & Lainey (2007). It can be shown that in the particular case of a = —1 
the parameter S coincides with At. In realistic situations, a differs from —1, while the 
parameter S is related to the time lag in a more sophisticated way (Ibid.). 

To show how (72 - 74) stem from the above Assertions, keep for the time being a — —1 . 
Also recall that the torque is despinning (so 9 > n), and that for the averages over time 

(^■) = ^ , (76) 

C = ^ Mpianet Fl? being the maximal moment of inertia of the planet. (For a homogeneous 
spherical planet, ^ = 2/5.) Then plug (75) into (59) and average the torque:^^ 

2 6 

^ _ iGM,^ he ( _ .J fl_ ^ +0{iyQ)+0(Q-')+0(en/Q\) = 
- i^^LM e{^) + i^il^tM ( , :J ) + oayq) + 0(Q-^) + 0(en/Q\) (77a) 



^ -(1-e^) - (1 + ecos.) d. 



R 



R Cb ZTT Jo 

where the absolute error 0{en/Q'^x) emerges due to an uncertainty in the definition of the 
overall quality factor Q employed in MacDonald's model. 

Evaluation of the above integrals is trivial and indeed leads to (72 - 74), the constant being 

^ ^ SCmI k2S R ^ 3n''M'^,k2At r' ^ Sum'^, k, r' n 

C R a6 e Mpianet {Mpi^net + M^at) Mpi^net {Mpianet + M^at) X 

where we used the fact that for a = —1 the rheological parameter S is simply the lag At. 

It should also be added that, since (77b) contains a relative error 0(Q~^) , the usefulness 
of the and terms in (73 - 74) depends on the values of the eccentricity and the quality 

As explained in the paragraph preceding formula (59), substitution of sine with 1/Q in the expression 
for torque generates a relative error 0{Q~'^) , i.e., an absolute error 0{Q~^) . Instead of inserting (75) into 
(59), one may directly use (53). Still, approximation of sine with e will entail, in (77) and its coroUaria, a 
relative error 0{Q~'^) and an absolute error 0{Q^^) . The situation will become more complicated in the 
special case of low values of the quality factor. See the Appendix below 
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factor. If, for example, Q = 70 , then the terms become unimportant for e < 0.12 , while 
the terms become unimportant for e < 0.24 . 

To draw to a close, we would mention that besides the above formula (72), in the literature 
hitherto we saw its sibling, an expression derived in a similar way, but with Assertion II rejected 
in favour of treating Q as a frequency- independent constant. The result of this treatment 
suffers an incurable birth trauma - the incompatibility between the frequency-independence 
of At and the frequency-independence of Q . 



10.2 Calculation based on the Darwin torque 

The following alternative derivation is based on the same Assertions (I - III) and, naturally, 
leads to the same results. The idea is to calculate the despinning rate not in terms of the 
MacDonald torque, but in terms of the Darwin torque, keeping the eccentricity-caused relative 
error at the level of O(e^) . 

To keep the inclination-caused relative error at the level of O(i^) , we still assume, in (34), 
that 1 = 2, m = 2 , p = . As for the the values of q , wc keep only the ones giving us 
terms of order up to e"^ , inclusively. Besides, we assume the phase lags to be small, so that 
^^^^impq = ^impq + ^i^^) = ^impq + 0{Q~^) . Under all these presumptions, the constant 
part of the tidal torque can be approximated with 

3 ^ 2 

^1=2 = oGMl.R'a-'k, G sine^^^^ + 0{e'/Q) + 0{e/Q) (79a) 

q=-2 

= ^GMl,R^a-'k, G,oq'22oq + ^^e'/Q) + O(zVg) -h 0{Q-'^) , (79b) 

q=-2 

where, according to the tables (Kaula 1966), 

= , = ^ _ ^ + 0(e6) , = 1 - 5e2 + - + 0(6^) , 

20-2 ' 20-1 4 16 200 8 ' V ; > 

^ 49 861 4 6 =—e^ + 0{e') , 

201 4 16 ' \ I ■> 202 4 ' V y > 

and, according to formula (24), 

e = [-2 9) At , e = [-29 + n) At , e =(-2^ + 2n)At 

220 -2 ^ ' 220 -2 ' 220 -1 ^ ' 220 -1 ' 2200 ^ ' 2200 



2 ^ + 3 n) At^ 



2^ + 4n) At 



2201 ^ ' 2201 ' 2202 ^ ' 2202 

Provided the quality factor scales as inverse frequency, all the time lags are the same constant 
At , so the above formulae all together entail, in the case of nonresonant prograde spin: 



^(l + ye^ + ^e^ + 0(e^) 



^2) 



27 2 573 



+ o(^Vg) + o(g-3) , 



which coincides with (72 - 74) to the order , inclusively, provided we substitute X2200 iiistead 
of X in the expression (78) for /C . 
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10.3 Rheologies different from Q ~ 1/x 

A part and parcel of both afore-presented methods was the assertion of all the time lags Ai^^^^ 
being equal. In reality, the time lags vary from one harmonic to another. 

Any particular functional form of the dependence At{x) fixes the rheology: for example, 
the frequency-independence of At constrains the value of the exponential a to —1 (while 
the parameter S becomes simply At ) . However, for an arbitrary a 7^ — 1 the lags will read 
(Efroimsky & Lainey 2007): 

-a -(a+1) 

At, ^ 8 X 83 

While the MacDonald approach cannot be generalised to a ^ — 1 , the Darwin-Kaula- 
Goldreich method can be well combined with (83). To this end, we shall insert (80 - 81) 
and (83) into (79a), and shall also employ the evident formulae 



a a 



cote, I Q, Ex, 

Impq ' Impq Impq /0/i\ 

cose = , = , = = , , (84) 

1 + 'impq Jl + g, Jl + S X, 



Impq y Impq 



sgntJ, sgncj, sgncc;, 

Irnpq Impq Impq 



sme, = sm e, sgno;, = — 1= , , 

" " Vi + ^^^Q^ ,/i + /" 

y Impq y Impq 

with uj, given by (25). and |e, I assumed (for reasons explained in the Appendix) not to 

Impq -J \ y . I luipq ' ^ ± i. / 

approach too close to 7r/2 . This will give us the following expression for (the constant part 
of) the deceleration rate of a non-resonant prograde spin: 



.. ^ T ^ 3 G Ml, k. 



C 2 a? ^MpianetR^ 



e ^ S \2e - n\ 

— sgn(2 6' — n) 



I 2 63 A , ■ , |2^ - 2nr 

+ 1 - 56^ + — sgn(2 9 - 2n) 



8 V ~ ' 1 + e^'' \2Q - 2n''" 



,'49 2 861 A , • , E^" \2Q - 3nr 
+ — - — eM sgn(2^ - 3n) 



4 16 



1 + 12^ - 3n|'" 



289 4 • , , \2e - 4nr 

+ — sgn(2 6' - 4n) 



1 + £: 26* - 4n 



+ O(zVQ) + O(eVQ) 



Be mindful, that a naive substitution of the formula (85) for sin e^^^^ into (79a) would result in 
an expression for the torque, attaining its maxima on approach to resonances (for a positive a), 
an evidently unphysical behaviour. As explained in section 9, there exists a profound physical 
reason, for which the actual multiplier in (79a) must be not sine^^^^ but: sine^^^^ '^'^^^impq ■ 
Mathematically, the presence of the cosine is irrelevant unless X^^^^ and Q^^^^ approach zero. 
If however X^^^^ becomes very small (i.e., if we approach a resonance), it is this long-omitted 
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(though known yet to Darwin 1879) cosine multipher that saves us from the unphysical maxima 

- see section 9 above. 

Under the extra assumptions'^^ of le^^^gl <^ 1 and n <ti 9 , formula (86) simphfies to 



-^(l + ye^ + ^e^ + 0(e- 



+ n ( 1 + ( y - 6q; I e- 



105 363 ,4 ^, 6^ 
— —a ) + 0(e6 



0{i'/Q)+0{Q-') + 0{n/e) 



IC 



n ( 1 + ( y - 6a ) 



where the overall factor /C is given by 



R 



3 

R n 



C Mpianet {Mpianet + M^at) C ^2200 {Mplanet + M^at) X 



'-2200 



(89) 



an expression identical to (78), except that it contains Af22oo > Q2200 ' X2200 instead of 
At , Q , and x i correspondingly. 

Were a equal to — 1 , sum (87) would coincide with (82), provided 9 > 2n (but not 
otherwise!). For realistic mantles and crusts, though, the values of a will, as pointed above, 
reside within the interval 0.2 — 0.4 (closer to 0.2 for partial melts). 



The smallness of enables one to employ (79b) instead of (79a). Then the multipliers 



s n " 

in (86) become e 



Impq 

J2a 2a 



Impq 



At, 



Impq 2200 ~Kt 



Impq 



2200 



^Impq ^^2200 ^i^pq 1 X 



2200 
Impq 



a+1 



^2200 ^*2200 ^Impq \ X 



2200 
Impq 



>^2200 ^*2200 «g°^^„^p, ' ^ + 



^Impq 



^2200 



X. 



2200 



. Specifically, 



S2og =«S'i^' 



At 



220g ""2200 '^'^2200 



^ , ^220g -•^2200 
^2200 



sgn w, 



220(7 2200 ^2200 



Xo 



_ ^ -^220^ "^2200 



'-2200 



^§'^'^2209^*2200 + 



2200 ^ ^22oq ' latter approximation being legitimate only under the 



condition of X22Qq ^ ^2200 ^ -'^2200 ' ^^^"^^ turns out to be equivalent to n ^ ^ . For example, 



^^2209 ~ X2200 ^ (-26l+4n)-(-26>+2n) ^ 



2n 



^^2200 



-2 6* +2n 



2 6* +2n 



So approximating £220(7 ' 



9 = — 2,— 1,0,1,2, we arrive at formula (87). 

There exists one more reason to keep n much smaller than 9 in (86 - 89). We derived (86 - 89) by inserting 
the customary relation (36) into (84 - 85). As explained in the Appendix below, (36) becomes invalid near 
spin-orbit commensurabilities. Indeed, at each commensurability a certain tidal harmonic becomes nil - see 
formula (26). According to (60), the appropriate Q, too, becomes nil. In this situation, one has to rely not 
on (36) but on a more general formula (105). The latter formula however entails vanishing of the appropriate 
component of the tidal torque on crossing the commensurability - see (113 - 114). This is why in (87 - 89) and 
even earlier, in (86), we should stay away from the commensurabilities 6 = n/2 , 9 = n , 9 = 3n/2 , or 9 = 2n . 
So we better keep n <^ 9 . 
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11 Conclusions 



In the article thus far we have provided a detailed review of a narrow range of topics. Our 
goal was to punctiliously spell out the assumptions that often remain implicit, and to bring to 
light those steps in calculations, which are often omitted as "self-evident" . 

This has helped us to demonstrate that MacDonald-stylc formula (57) for the tidal torque 
is valid only in the zeroth order of en/ Qx , while its time-average is valid only in the first order. 
These restrictions mean that in the popular expressions for tidal despinning rate the terms with 
higher powers of e become significant only for large eccentricities. Their significance is limited 
even further by the error 0{Q~^) emerging when the sine of the phase lag gets approximated 
by the inverse quality factor - see formula (59) and the paragraph preceding it. 

We have demonstrated that in the case, when the inclinations are small and the phase 
lags of the tidal harmonics are proportional to the frequency, the Darwin-Kaula expansion is 
equivalent to a corrected version of the MacDonald formalism. The latter method rests on the 
assumption of existence of one total double bulge. The necessary correction to MacDonald's 
approach would be to assert (following Singer 1968) that the phase lag of this integral bulge is 
not constant, but is proportional to the instantaneous synodal frequency 2(i> — 6*) , where u and 
9 are the true anomaly and the sidereal angle. Any rheology different from this one will violate 
the equivalence of the Darwin-Kaula and MacDonald descriptions. It remains unexplored if 
their equivalence is violated also by setting the inclination high. 

We have demonstrated that no "paradoxes" ensue from the frequency-dependence Q ~ > 
with q; = 0.3 ± 0.1 , found for the mantle. 

We have investigated the limitations of the popular formula interconnecting the quality 
factor Q and the phase lag e. It turns out that for low quahty factors (less than 10), the 
customary formula Q = cot |e| should be substituted with a far more complicated relation. 

Finally, we examined two derivations of the popular expressions (72 - 74), and have pointed 
out that these expressions have limitations related to the frequency-dependence of the quality 
factor. First, dependent upon the values of e and Q, the high-order terms in these expressions 
may become significant only for large eccentricities. Second, the expansion of the deceleration 
rate in even powers of e will be different if At is frequency-dependent (which is the case for 
solid materials). These two circumstances do not necessarily disprove any major result achieved 
in the bodily-tide theory. However, some coefficients may now have to be reconsidered. 

For the reahstic rheology of terrestrial bodies, the despinning rate, in the absence of tidal 
locking, is given by our formulae (86 - 89). 
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Appendix. 



The lag and the quaUty factor: is the formula Q = cot |e| universal? 

The interrelation between the quahty factor Q and the phase lag e is long-known to be 

Q = cot |e| , (90) 

and its derivation can be found in many books. In Appendix A2 of Efroimsky & Lainey(2007), 
that derivation is reproduced, with several details that are normally omitted in the literature. 
Among other things, we pointed out that the interrelation has exactly the form (90) only in 
the hmit of small lags. For large phase lags, the form of this relation will change considerably. 

Since in section 9 of the current paper we address the case of large lags, it would be worth 
reconsidering the derivation presented in Efroimsky & Lainey (2007), and correcting a subtle 
omission made there. Before writing formulae, let us recall that, at each frequency x the 
spectrum of the deformation, the quality factor (divided by 2 tt ) is defined as the peak energy 
stored in the system divided by the energy damped over a cycle of fiexure: 

where AE(.yde{x) < as we are talking about energy losses. 

An attempt to consider large lags (all the way up to |e| = 7i/2) sets the values of Q/27i 
below unity. As the dissipated energy cannot exceed the energy stored in a free oscillator, the 
question becomes whether the values of Q/27i can be that small. To understand that they 
can, recall that in this situation we are considering an oscillator, which is not free but is driven 
(and is overdamped). The quality factor being much less than unity simply implies that the 
eigenfrequencies get damped away during less than one oscillation. Nonetheless, motion goes 
on due to the driving force. 

Now let us switch to the specific context of tides. To begin with, let us recall that the 
dissipation rate in a tidally distorted primary is well approximated by the work that the 
secondary carries out to deform the primary: 

E ^ - [ p V ■ VW (92) 



p , V , and W denoting the density, velocity, and tidal potential inside the primary. The 
expression on the right-hand side can be transformed by means of the formula 

p V -vw = V -ipyw) - w V -Vp - wv -ip v) = v-ipvw)-wv-vp + w^ 

at 

where the W V -V p and dp/dt terms may be omitted under the assumption that the primary is 
homogeneous and incompressible. In this approximation, the attenuation rate becomes simply 

E ^ - j V ■ {p V W)d^x ^ - j pW V ■ n dA , (94) 

n being the outward normal to the surface of the primary, and dA being an element of 
the surface area. It is now clear that, under the said assertions, it is sufficient to take into 



We are considering flexure in the linear approximation. Thus at each frequency x the appropriate energy 
loss over a cycle, A£Jcj/cie(x) > depends solely on the maximal energy stored at that same frequency, Epeakix) ■ 
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account only the radial elevation rate, not the horizontal distortion. This way, formula (92), 
in application to a unit mass, will get simplified to 

C standing for the vertical displacement (which is, of course, delayed in time, compared to 
W ) . The amount of energy dissipated over a time interval {to , t) is then 

We shall consider the simple case of an equatorial moon on a circular orbit. At each point 
of the planet, the variable part of the tidal potential produced by this moon will read 

W ^ Wo cosxt , (97) 

the tidal frequency being given by 

X = 2\n - ujp\ . (98) 

Let g denote the surface free-fall acceleration. An element of the planet's surface lying beneath 
the satellite's trajectory will then experience a vertical elevation of 

W 

C = h2— cos(xt - |e|) , (99) 

/i2 being the corresponding Love number^^, and |e| being the positive^^ phase lag, which for 
the principal tidal frequency is simply the double geometric angle S subtended at the primary's 
centre between the directions to the secondary and to the main bulge: 

|e| = 25 . (100) 

Accordingly, the vertical velocity of this element of the planet's surface will amount to 



u — ( — —h2X~~ sin(xi — |e|) = ~ ^2 X—~ (sin^i cos|e| — cosxt sin|e|) . (101) 



The expression for the velocity has such a simple form because in this case the instantaneous 
frequency x is constant. The satellite generates two bulges (on the facing and opposite sides of 
the planet) so each point of the surface is uplifted twice through a cycle. This entails the factor 
of two in the expression (98) for the frequency. The phase in (100), too, is doubled, though the 
necessity of this is less evident, - see footnote 4 in Appendix to Efroimsky & Lainey (2007). 



For a homogeneous incompressible body, k2 = (3/5)/i2 , for which reason (99) and the subsequent equations 
with /i2 can equally be written as proportional to ■ The formulation employing ^2 is more fundamental, 
as it can, in principle, be generalised to a compressible body of a radially-changing density. Indeed, whatever 
the properties of the primary arc, the dissipation rate in it is equal to the rate of change of the primary's spin 
energy plus the rate of change of the orbital energy. Both the latter and the former are proportional to ^2 . 

Were we not considering the simple case of a circular orbit, then, rigorously speaking, the expression for 
W would read not as Wo cosxt but as Wo cosw^-^^jt, the tidal frequency taking both positive and 

negative values, and the physical frequency of flexure being x = \'^udai I ■ Accordingly, the expression for C, 
would contain not cos(xi — |e|) but cos(a;j-jj^;t — e) • As we saw in equation (24), the sign of e is always the 
same as that of . For this reason, one may simply deal with the physical frequency x = \'^tidai I ^^'^ 

with the absolute value of the phase lag, |e| . 
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The energy dissipated over a time cycle T — 27r/x, per unit mass, will, in neglect of 
horizontal displacements, be 



AE 



cycle 



[■(- 



dW\ , / , Wo\ dWo f'=^ , . , , . , ,^ , 

— — \ at = — \ — h2 X — — / cosvt (sm yt cos e — cosvt sm e j at 
dr J V g J dr Jt^o 



Xt = 27T 



= - ^2 X — sin |e| - / cos^x^ d{xt) = - /12 — tt sm |e| 

g dr X Jxt = o g dr 

while the peak energy stored in the system during the cycle will read: 



(102) 



E 



peak 



T/4 



U 



dt = 



/i2 X — — / cosxt (smxt cos e — cosx^ sm e j at 
g ) dr 



g dr 



cos e 



X 



cosxt sinx^ d{xt) 



xt=\A 



sm e 



X 



cos^X^ d{xt) 



xt = \A 



In the appropriate expression in Appendix Al to Efroimsky & Lainey (2007), the lower limit 
of integration was erroneously set to be zero. To understand that in reality integration over xt 
should begin from |e|, one should superimpose the plots of the two functions involved, cosx^ 
and sin(xt — |e|). The maximal energy gets stored in the system after integration through the 
entire interval over which both functions have the same sign. Hence x^ — kl ^s the lower limit. 
Evaluation of the integrals entails: 



E. 



peak 



whence 



Wo dWo 



g dr 



cos e 
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1 cos 



|e| — 2 ^ ■J ~ kl ) sin |e| ~ ( I' ~ 1^1 ) kl 



As can be seen from (105), both the product sine cose and the appropriate component of the 
torque attain their maxima when Q — 1 — tt/A . 

Usually, |e| is small, and we arrive at the customary expression 

Q-^ = tan|e| + O(e^) . (106) 

In the opposite case, when Q — > and |e| 7r/2, it is convenient to employ the small difference 

|e| , (107) 



^ 2 



in terms whereof the inverse quality factor will read: 
cot C 1 1 



^ cot^ tsm^ 



arctan z 



1 y3 

3 ^ 



1 - ^ + 0(z^ 



, (108) 



. (103) 



(105) 
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where z = tan,^ and, accordingly, ^ = arctanz = z — ^ z^ + ^ z^ + 0{z'^) . Formula 
(108) may, of course, be rewritten as 



1 - ^z' + 
5 



(109) 



or, the same, as 



z = (3g)V3 



1 + ^z' + 0{z') 
5 



(110) 



While the zeroth approximation is simply z — {SQY^^ + 0{Q) , the first iteration gives: 

1 



tan^ = z — (Sg) 



1/3 



1 + I (3g)^/3 + 0{Q'/') 
5 



, (111) 



with q = (Sg)^/^ playing the role of a small parameter. 

We now sec that the customary relation (106) should be substituted, for large lags, i.e., for 
smalP^ values of Q , with: 



tan|e| = (Sg)"^/^ 



1 - I m'^' + 0{Q'/') 
5 



(112) 



The formula for the tidal torque contains a multiplier sin e cos e , whose absolute value can, 
for our purposes, be written down as 



sin|e| cos|e| =cos^ sin^ 



whence 



tan,^ 



i + 5g' + 0(g^ 



l + tan2{ l + g2[i + 0(g2 



^ = (3g)V3 



sme cose 



±(3g) 



1/3 



l_|(3g)2/3 + 0(g4/3) 

5 



l_^(3g)2/3 + 0(g4/3) 

5 



(114) 



an expression vanishing for g — >^ . Be mindful that both e22oo ^^^1 the appropriate compo- 
nent of the torque change their sign on the satellite crossing the synchronous orbit. 
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